Quasiparticle Resonant States Induced by a Unitary Impurity in a d-Wave 

Superconductor 
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The quasiparticle resonant states around a single nonmagnetic impurity with unitary scattering in 
a d-wave superconductor is studied by solving the Bogoliubov-de Gennes equations based on a t-J 
model. Both the spatial variation of the order parameter and the local density of states (LDOS) 
around the impurity have been investigated. We find: (i) A particle-hole symmetric system has 
a single symmetric zero-energy peak in the LDOS regardless of the size of the superconducting 
coherence length £o; (h) For the particle-hole asymmetric case, an asymmetric splitting of the zero- 
energy peak is intrinsic to a system with a small value of &;f£o- 



PACS numbers: 74.25. Jb, 74.50.+r, 73.20.Hb 

It is now well established [Q that high-T c supercon- 
ductors (HTSC's) have essentially a d x 2_ y 2-wave pair- 
ing symmetry. In conventional s-wavc superconductors, 
nonmagnetic impurities affect neither the transition tem- 
perature nor the superfluid density as dictated by the 
Anderson theorem 0. But in a d-wave superconductor 
(DWSC) with nodes of the energy gap, such impurities 
can cause a strong pair-breaking effect || . Recently, the 
local electronic properties in the immediate vicinity of an 
isolated non-magnetic impurity in a DWSC has become 
the topic of increased investigation |4|-[l(| , as these prop- 
erties may provide a distinctive signature for the pairing 
symmetry. It has been theoretically predicted by Bal- 
atsky, Salkola and co-workers J?],|| that, in a DWSC, a 
single nonmagnetic impurity can generate quasiparticle 
resonant states at subgap energies. They showed that, for 
a moderately strong impurity, an asymmetry of the reso- 
nance peak near the Fermi energy is induced by the fact 
that the impurity locally breaks the particle-hole sym- 
metry. However, their theory says that increasing the 
impurity strength pushes the resonance peak toward the 
Fermi level, so that, in the unitary limit, the resonance 
occurs right on the Fermi level, and only a single symmet- 
ric zero-energy peak (ZEP) occurs in the LDOS near the 
impurity. It has also been shown by a finite-size diago- 
nalization |]l0f that, in the unitary limit, the lowest eigen- 
values are essentially zero, indicative of the appearance 
of zero-energy states (ZES's). Note that, in Ref. [[To|| , 
the chemical potential /z was taken to be at the center 
of the tight-binding energy band (i.e., fi = 0), so that 
the system has a particle-hole symmetry. This symme- 
try is also upheld in the continuum-theory treatment of 
impurities where the self-consistent i-matrix approx- 
imation is employed. A question which arises naturally 
is whether, in the unitary limit, the "ZEP" in the LDOS 
due to the "ZES's" has an asymmetric splitting or not, 
when particle-hole symmetry is broken in the system. 
Recently, Tanaka et al. p5| concluded with their numer- 



ical study that such a splitting is still present, whereas 
Tsuchiura et al. |l6| made an opposite conclusion in their 
numerical study, and asserted that the system studied by 
Tanaka et al. was too small for their results to be reli- 
able. Experimentally, an asymmetric splitting is clearly 
observed by Yazdani et al. , whereas Hudson et al. Jp| 
observed only an off-centered peak with no indication of 
a splitting. Thus it appears important to settle the issue 
of whether a unitary non-magnetic impurity in a pure 
DWSC can indeed give rise to such an asymmetric split- 
ting in the "ZEP" , as it will decide whether experimental 
observation of this feature in HTSC's necessarily implies 
that these SC's do not have pure d-wave symmetry, or 
that the impurity is not in the unitary limit (in which 
case the asymmetry is tied to the sign of the impurity 
potential, which may well be a misleading conclusion). 

Based on a t-J model, this paper presents an extensive 
study on the electronic states around a unitary single-site 
impurity in a DWSC. The spatial variation of the su- 
perconducting order parameter (OP) near the impurity, 
including an induced s-wave component, is determined 
self-consistently. By investigating the sensitivity of the 
LDOS on both fi and £o, we nn d: (i) When fi — 0, so 
that the system is particle-hole symmetric, a single ZEP 
occurs in the LDOS spectrum which is symmetric with 
respect to zero energy, regardless of the size of £o; (ii) As 
the particle- hole symmetry is broken by letting fi =/= 0, 
a critical value j c exists, which is larger for larger 
so that for 7 = kp£,o < 7c the "ZEP" exhibits an asym- 
metric splitting. (Here Uf is the Fermi wavevector.) [T^ ] 
Thus we find that for a particle-hole asymmetric sys- 
tem, a sufficiently small coherence length can cause the 
"ZEP" to exhibit an asymmetric splitting. Treating such 
a system by the self-consistent i-matrix approximation, 
which restores the particle-hole symmetry, will then lose 
this feature and be misleading in this respect. 

We consider a t-J model Hamiltonian defined on a two- 
dimensional square lattice: 
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<ij> <ij> 

where the Hilbert space is made of empty and singly- 
occupied sites only; summing over (ij) means summing 
over nearest-neighbor sites; rii = ^ CT C; CT Ci CT is the elec- 
tron number operator on site i; Si is the spin-i operator 
on site i; and J > gives the antiferromagnetic superex- 
change interaction. As in Ref . yj| , we have also included 
a direct nearest-neighbor interaction term. W — and 
J/4 correspond to two versions of the standard t-J model. 
This term is introduced to adjust the magnitude of the 
resultant <i-wave OP. The scattering potential from the 
single-site impurity is modeled by U\ — Uq8u with / the 
index for the impurity site. The slave-boson method |HJ 
is employed to write the electron operator as c\ a = 
where j\ a and b\ are the operators for a spinon (a neutral 
spin-i fermion) and a holon (a spinless charged boson). 
Due to the holon Bose condensation at low temperatures, 
the quasiparticles are determined by the spinon degree of 
freedom only. Within the mean-field approximation, the 
Bogoliubov-de Gennes (BdG) equations are derived to be 

?(3 £)(?)-*(?)■ <» 

with 

ffy =-[t6+(t + W) Xii ]S. +s . + (Ui - M )<5y . (3) 

Here u™ and i?" are the Bogoliubov amplitudes corre- 
sponding to the eigenvalue E n ; 5 and Xy are the doping 
rate and the bond OP, respectively; and 5 are the unit 
vectors ±x, ±y. The resonant-valence-bond (RVB) OP 
Ay, xy, and 8 are determined self-consistently: 

Ay = Y>>T + u ? v ?*} t^h(E n /2k B T)6 i+s . , 

n 

(4) 

xy = £kx/(£«) + v >r^ - f( E n)}} > (5) 

a 

and 

S=1 ~W a EiKI 2 /^™) + K! 2 [l - f(En)}} , (6) 
i,n 

where kg is the Boltzmann constant; f(E) = 
[exp(E/kBT) + is the Fermi distribution function; 
and N a = N x x N y is the number of lattice sites. The 
BdG equations is solved fully self-consistently for the 
bulk state first. We then fix the values of 6 and x an d 



solve the BdG equations in the presence of a single impu- 
rity with the self-consistent <i-wave RVB OP. The ther- 
mally broadened local density of states (LDOS) is then 
evaluated according to: 

Pi (E) - -2]T[|<| 2 /'(£„ -E) + \vVff(E n + E)] , 

n 

(7) 

where a factor 2 arises from the spin sum, and f'(E) = 
df{E)/dE. The LDOS pi(E) is proportional to the lo- 
cal differential tunneling conductance which can be mea- 
sured in a scanning tunneling microscope/spectroscopy 
(STM/S) experiment ||. 

In the numerical calculation, we construct a super- 
lattice with the square lattice N x x N y as a unit su- 
percell. As detailed in Ref. |2C|] , this method can pro- 
vide the required energy resolution for the possible res- 
onant states. Throughout the work, we take the size 
of the unit supercell N a = 35 x 35, the number of 
supercells iV c = 6x6, the temperature T = 0.01J, 
and the single impurity potential in the unitary limit 
Ua = 100J. The values of the other parameters — /z, 
W, and t, are varied in order to investigate the elec- 
tronic states around a single impurity for various ways 
to bring about particle-hole asymmetry. The obtained 
spatial variation of the d- wave and the induced extended- 
s-wave OP components around the impurity, which are 
defined as A d (i) = \[A x (i) + A_ 4 (i) - A y (i) - A_g(i)], 
and A s (i) = \ [A* (i) + A_ £ (i) + Ag (i) + A_ fl (i)] , is sim- 
ilar to Fig. 1 of Ref. ||. These OP component's have 
the following characteristics: The d-wave component de- 
creases continuously to zero from its bulk value as the 
impurity site is approached, in the scale of the coher- 
ence length £o = hvF /ft A max , with the depleted region 
extending farther in the nodal directions if £o is larger. 
Here A max = 4Ao with Ao the bulk value of the d-wave 
OP defined in the real space on a nearest neighbor bond, 
and vf is the Fermi velocity. The s-wave component is 
zero at the impurity site and also at infinity. It has line- 
nodes along the {110} and {110} directions, and changes 
sign across any nodal line. Unlike the pairing state at 
a {110} surface of a DWSC, which can break the time- 
reversal symmetry, the pairing state near a single impu- 
rity conserves time-reversal symmetry. This difference 
can be understood from the Ginzburg-Landau (GL) the- 
ory f22| , in that a mixed gradient term favors the d- and 
induced s-wave OP components to be in phase, but it 
vanishes near a {110} surface, whence the fourth order 
s-d coupling term can establish an s + id pairing state. 

Figure |l| shows the LDOS as a function of energy on 
sites one and two lattice spacings along the (100) direc- 
tion from the impurity and on the corner site of the unit 
supercell. The values of the parameters are labeled on 
the figure. Note that the LDOS at the corner site has 
recovered the bulk DOS, by exhibiting a gaplike feature 



2 



with the gap edges at ±A max . This resemblance indi- 
cates that the unit cell size and the number of unit cells 
are large enough for uncovering the physics intrinsic to 
an isolated impurity. As shown in Fig. El we find that 
the LDOS spectrum near the impurity is highly sensitive 
to the position of /i within the energy band, and the pa- 
rameter 7. In Fig. 0(a), fx = ||| and 7 = 0.80, a single 
ZEP occurs in the LDOS on the nearest-neighbor site of 
the impurity, similar to the prediction of the continuum 
theory [0,|| and the eigenvalue calculation in Ref. [ pj . 
In addition, as a reflection of the particle-hole symmetry, 
the whole LDOS spectrum is symmetric about E = 0. 
We have also studied the cases (not shown) with the 
same \x = and t — 4J but with W = and W — 0.5J 
(corresponding to 7 = 0.27 and 2.0), and found that 
the above feature remains unchanged, which allows us 
to conclude that as long as the system is particle-hole 
symmetric, only a single symmetric ZEP exists for all 7. 
When /1 is not zero, the system is particle-hole asym- 
metric, and the LDOS spectrum becomes asymmetric. 
[See Fig. 0(b)-(g).] In Fig. 0(b)-(e), n = -0.32 J is fixed, 
and W and t are varied in order to change 7. For a 
large 7 = 91.7, we see a single ZEP in the LDOS. [See 
Fig. |(b).] When 7 is lowered to 16.5, the "ZEP" be- 
gins to evolve into a double-peaked structure with the 
E > peak having the dominant spectral weight over 
the E < peak. For a further decreased 7 = 6.05, the 
spectral weight of the peak at E < is enhanced. (See 
Fig. 0(d).) As seen in Fig. 0(e), this enhancement be- 
comes even more pronounced when [i is made close to 
the edge of a very narrow energy band so that 7 becomes 
as small as 2.85. When /j, = —0.16 J, we only observe 
a single ZEP although 7 is as small as 5.7 [for Fig. 0(f)] 
and 2.9 [for Fig. 0(g)], except that a tendency of the split- 
ting can be identified in the latter case. This tendency of 
the splitting has been observed clearly in STM tunneling 
spectroscopy measurements (see Fig. 4(A) of Ref. [^2|). 
It should be emphasized that the ZEP splitting obtained 
here has a different origin from that found by Tanaka et 
al. ]lq| . We have re-examined their results by choosing 
the same parameter values and the system size (18 x 18). 
When the LDOS spectrum is displayed in a wide energy 
landscape, many split DOS peaks appear with no well- 
defined gaplike feature identifiable. But as the system 
size is enlarged by the supercell technique, the calculation 
only shows a single ZEP in the LDOS, which indicates 
that the splitting of ZEP obtained in Ref. pa] is indeed 
due to the size effect. On the other hand, we have also 
calculated the excess charge distribution due to the pres- 
ence of the impurity (oc Sni = (rii) — hq, where no is the 
average particle occupation on each site for the bulk sys- 
tem). We find that this distribution is anisotropic, with 
its magnitude having tails extending along the nodal di- 
rections [See Fig. |^j. Because Fig. is obtained with the 
parameter values given in Fig. 0(e) which gives a small 
7 (= 2.85) value, the exhibited tail is short. A similar 



calculation with the model parameters given in Fig. Iff) 
(not shown) shows that the charge distribution is sim- 
ilar to that displayed in Fig. || except for a longer tail 
along the nodal directions due to the larger 7 (= 5.7). 
This similarity in the charge distributions for a split and 
a unsplit ZEP's disproves the assertion made in Ref. pl| 
that the local charge-density oscillation is the cause of the 
ZEP splitting. We mention in passing that we have also 
found that the excess charge density decays exponentially 
along the nodal directions instead of the r _2 -dependcncc 
from the impurity. But we do not think that this finding 
invalidates the assertion in Ref. ffi] that the wavefunction 
of the impurity resonant state has a 1/r decay along the 
nodal directions, which can lead to a long range inter- 
action between the impurities. However, we do believe 
that since we have obtained essentially the bulk density 
of states in several neighboring points near the corner 
of the supercell, the interaction between the neighboring 
impurities should be negligible in the cell size we have 
chosen to work with. Thus we believe that it is very 
unlikely that the splitting of the ZEP we obtain is due 
to this interaction. Since the s-wave OP component in- 
duced near the impurity is in phase with the dominant 
d-wave component, the splitting of the ZEP we found 
is not due to a local broken time-reversal symmetry. Fi- 
nally, as shown in Fig. 0(e), the splitting is also exhibited 
in a non-self-consistent calculation with a spatially uni- 
form bulk d-wave OP, showing that the suppression of 
the d-wave OP component, and the induction of the s- 
wave component, have little to do with the splitting. All 
of these points lead us to the conclusion that, for the 
particle-hole asymmetric case, the splitting of the ZEP is 
intrinsic to the system with a short coherence length, and 
the critical value j c , below which the ZEP is split into 
an asymmetric double-peak, is simply a reflection that 
the system has reached a critical extent in its deviation 
from particle-hole symmetry. We thus propose to under- 
stand these results qualitatively as follows: The "ZES's" 
induced by a unitary non-magnetic impurity have essen- 
tially the same physical origin as the "midgap states" pre- 
dicted to exist on the surfaces/interfaces of a DWSC [§3 . 
Their existence is implied topologically by the Atiyah- 
Singer index theorem |25[| , which applies to particle-hole- 
symmetric Dirac-like operators. When this symmetry 
is mildly broken the midgap states are expected to still 
exist but no longer exactly "midgap" . The BdG equa- 
tions become Dirac-like equations only under the WKBJ 
approximation (which is a part of the self-consistent t- 
matrix approximation), the error of which is measured 
by the parameters and 7 . For fi = 0, the system 
has exact particle-hole symmetry for all 7. Thus smaller 
\fi\ should imply smaller 7 needed to reach the same de- 
viation from particle-hole symmetry, hence a smaller j c 
below which an asymmetric splitting of the ZEP appears. 

In summary, we have presented an extensive study on 
the quasiparticle resonant states induced by a unitary 
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non-magnetic impurity in a DWSC. The results have clar- 
ified some conflicting conclusions in the literature, and 
should be of value for the proper analysis of the STM/S 
results obtained on HTSC's around an isolated impurity. 
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FIG. 1. Local density of states as a function of energy on 
sites one (solid line) and two lattice (dashed line) spacings 
along the (100) direction away from the impurity, and on the 
corner site (short-dashed line) of the unit cell. The parameter 
values have been correspondingly labeled on each panel. Also 
shown in the (e) panel is the local density of states on the site 
nestest-neighbor to the impurity obtained with a pure bulk 
d-wave order parameter (dotted line). 



FIG. 2. Spatial variation of the charge distribution around 
the impurity with the parameter values given in Fig. ^(e). 
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FIG. 1 Zhu et al. 




FIG. 2 Zhu et al. 



